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ABSTRACT 

The recent discovery that faint gamma-ray bursts are stretched in time 
relative to bright ones has been interpreted as support for cosmological distances: 
faint bursts have their durations redshifted relative to bright ones. It was 
pointed out, however, that the relative time stretching can also be produced by 
an intrinsic correlation between duration and luminosity of gamma-ray bursts 
in a nearby, bounded distribution. While both models can explain the average 
amount of time stretching, we find a difference between them in the way the 
duration distribution of faint bursts deviates from that of bright ones, assuming 
the luminosity function of gamma-ray bursts is independent of distance. This 
allows us to distinguish between these two broad classes of model on the 
basis of the duration distributions of gamma-ray bursts, leading perhaps to 
an unambiguous determination of the distance scale of gamma-ray bursts. We 
apply our proposed test to the second BATSE catalog and conclude, with some 
caution, that the data favor a cosmological interpretation of the time dilation. 
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1. Introduction 

From the data obtained by BATSE on faint gamma-ray bursts, combined with data 
from earlier satellites (notably PVO), we know the characteristics of the radial and angular 
distribution of gamma-ray bursts quite well: we are at or very near the center of the 
distribution of known gamma-ray bursts, and the bursts occur at a uniform rate density 
out to some characteristic distance, after which they become rarer. The big unknown is this 
characteristic distance, and the two most popular distance models place it in the very outer 
regions of our Galaxy (Galactic corona models) or at a fair fraction of the distance to the 
horizon of our universe (cosmological models). In principle, there is a third, 'Ptolemean' 
possibility: all gamma-ray bursts can be at the same substantial distance in a shell, and the 
shape of their flux distribution is simply a reflection of the true luminosity function, which 
by coincidence has a slope of —3/2 at the bright end. All other observed properties and 
the correlations between them are likewise intrinsic in this case. Since such models are not 
falsiflable until direct distance measurements to individual bursts become possible, we shall 
not discuss them further. 

The recently announced time dilation of a factor 2 in faint BATSE bursts (Norris 
et al. 1994, Lestrade et al. 1993; but see also Band 1994) has been adduced as evidence 
for the cosmological distance scale because the required redshift for the time delay is 
consistent with the redshift deduced from simple no-evolution, standard-candle models 
for cosmological gamma-ray bursts (Fenimore et al. 1993, Mao & Paczyhski 1992, Piran 
1992, Dermer 1992, Wickramasinghe et al. 1993). However, it was quickly pointed out 
that the effect can also be caused by an intrinsic anti-correlation between luminosity and 
duration of gamma-ray bursts. Such a correlation can be caused by many things, and 
relativistic beaming in particular has been discussed in some detail (Brainerd 1994, Mao & 
Yi 1994, Yi & Mao 1994). Another obvious candidate for introducing such a correlation 
is requiring a roughly constant total energy output, E oc L At. Any such anti-correlation 
will make it possible to get a lengthening of the average burst duration with decreasing flux 
even in Galactic corona models. 

In this paper, we present a way to easily visualize the interplay of distance and 
luminosity distributions in determining the observed log Ar(> F) — log of gamma-ray 
bursts (Sect. 2.). We use this to show that there is a generic difference between the two 
above models in the way the duration distribution of faint bursts compares with that of 
bright ones. We then apply a test based on this difference to the burst summary data as 
given in the second BATSE catalog, and flnd that a cosmological interpretation agrees best 
with the data (Sect. 3.). 
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2. Interpreting the flux distribution 



Both the cosmological explanation of the time dilation and the alternative local 
explanations that have been suggested are based on models in which the luminosity function 
of gamma-ray bursts is independent of distance, so we shall restrict ourselves to this case in 
the present paper. Some local models, such as those in which gamma-ray bursts originate 
on neutron stars that were shot out of our Galaxy, explicitly assume that gamma-ray burst 
properties change with increasing distance to us (Li & Dermer 1992). In these models, the 
correlation of duration with flux could be explained as an intrinsic property of the bursts. 
Just as in the above-mentioned Ptolemean model, explicit measurement of the distance 
to individual bursters is required to test such models, so we cannot discuss them in the 
present paper. The basic distinction we seek to explore is between models in which the 
distribution is bounded due to volume effects close to the horizon scale cj Hq (henceforth 
called 'cosmological models') and those in which the boundedness is due to a genuine 
decrease in volume density of bursters beyond some characteristic distance -Rcore ^ cj Hq 
(henceforth called 'local models'). Note that the latter encompass an enormous range 
of distance scales, from hundreds of astronomical units ('Oort cloud') up to hundreds of 
megaparsecs, and we offer no way of choosing among those. 

In Fig. 1 we illustrate how the flux distribution of gamma-ray bursts is shaped by 
their density and their luminosity function in the case of a local model. The central panel 
shows the luminosity of gamma-ray bursts as a function of their flux. For deflniteness, we 
assume that the luminosity function is a power law extending from Li to L^-, so all 
bursts in the central panel are between the horizontal dashed lines. The diagonal lines are 
lines of constant distance, with distance increasing from bottom right to top left. Again for 
deflniteness, we take a density distribution of the form n{r) oc [1 + (r/i?core)"] ^• 

There are now two possibilities for making the flux distribution (bottom panel of Fig. 1) 
agree with the observed log7V(> F) — log (a detailed discussion is presented in Ulmer & 
Wijers 1994). The flrst is that gamma-ray bursts are effectively standard candles (i.e. if 
Li ~ ^2 or if /3 > 2.5 or /3 < 1). Then the flattening of log N[> F) - log F for weak bursts 
simply reflects the decreasing density beyond -Rcore- The range of luminosities visible at 
each flux is narrow. The second possibility is that Li <C L2 and 1 < /3 < 2.5, in which case 
the slope below the break can simply be the slope of the (cumulative) luminosity function. 
The slope changes again to a value determined by the density at F-^. In this case, the range 
of luminosities seen at each flux level is quite broad. 



3. A generic diflFerence in the duration distributions 



- 4 - 



The durations of gamma-ray bursts range over decades, both at high and low fluxes, so 
the factor 2 duration difference between faint and bright bursts is a relatively small effect. 
We therefore assume that most of the duration range of bursts is intrinsic, and seek to 
explain a small non-intrinsic effect on top of that. 

3.1. Duration distributions in local models 

In local models, the longer average duration of faint bursts is explained by assuming 
that more luminous bursts tend to have shorter durations. As long as we are on the —3/2 
part of the flux distribution, the distribution of luminosity (and therefore duration) is 
independent of flux (Fig. 1). Thus, the duration distribution is the same for all fluxes above 
F2. But as soon as we go below F2, the effect of decreasing density will remove luminous 
(i.e. short) bursts from the sample, while leaving the weak (long) bursts unaffected. This 
effect progresses with decreasing flux until we pass Fi. Below Fi all bursts we see are 
beyond -Rcore! if the density distribution is an exact power law, we have once again reached 
a regime where all fluxes are equivalent, so the duration distribution does not change any 
further with decreasing flux (see Fig. 1). In case the low-flux slope is set by geometry 
because the luminosity function is narrow, this would still be true. But now Fi is effectively 
only just below F2, implying that the duration distribution changes its character over a 
narrow range in flux just below the break in the slope of log7V(> F) — log and remains 
constant thereafter. (Note that this offers a method, within the limited context of this 
model, to distinguish narrow luminosity functions from wide ones. Also bear in mind that 
the change must be substantial even if it occurs over a narrow range in flux in order that 
the factor of 2 increase in average duration be reproduced.) 

The implication is that for non-cosmological interpretations of the observed time 
dilation, one quite generically expects a change in the average duration of bursts to be 
the result of removing short bursts from the sample. That is to say, the long end of the 
gamma-ray burst duration distribution remains unchanged with decreasing flux, but the 
short end moves up, making the distribution signiflcantly narrower at low fluxes to increase 
its mean. 

3.2. Duration distributions in cosmological models 

In the case of cosmological distributions there is no need for wide luminosity functions 
or intrinsic correlations. Gamma-ray bursts can be near-standard candles occurring at 
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a constant rate per unit comoving volume and time (Mao & Paczynski 1992), and their 
other properties are independent of distance as well. In this case, the burst durations 
change simply due to redshift stretching. This is clearly different from the duration change 
discussed above because now the durations of all bursts increase, so the entire distribution 
shifts to longer durations and the width of the distribution is also increased by 1 + This 
means that the distribution becomes wider with decreasing flux, and its long end shifts up 
just as well as the short end. The relation between average duration and flux is flxed by the 
cosmological model one adopts. 

3.3. Comparing the distributions in practice 

In order to test our idea we used the second BATSE catalog (BATSE Team 1994), 
which is publicly available. Since it only contains summary data, we need to carefully select 
a sample that is less subject to known biases. This means we want to use peak fluxes in 
short time bins and use only bursts with durations much longer than those bins in order to 
avoid incompleteness at short durations and effects of fluence triggering. We therefore used 
the peak fluxes of bursts on the 64 ms time scale as a measure of brightness, and excluded 
all bursts shorter than 3 s from our sample. Also, we used the duration measure T^q rather 
than Tgo (Fishman et al. 1994), because the latter uses the far wings of the burst, which 
are subject to large errors especially in weak bursts (and are likely to lead to systematic 
effects in the relative duration measurements of faint and bright bursts). These issues can 
be dealt with using the full data (as was done by Norris et al. 1994), but not using the 
catalog. This leaves us with 187 bursts for which all necessary information is available. 
Following Norris et al., we compare the durations of the brightest bursts (top 10%) with 
those of a faint group, for which we choose those with count rates ranking from 10% to 
20% from the bottom (the faintest, low signal-to-noise group is avoided). The resulting 
histograms with 19 bursts in each group are shown in Fig. 2. The time dilation between the 
two groups is evident: the difference in mean log(r5o) is 0.34 it 0.11, signiflcant at the 3a 
level and consistent with the factor of 2 dilation found by Norris et al.. However, this was 
already known, and the issue here is which kind of time dilation is preferred. To address 
this, we note that the standard deviations of the distributions of log(r5o) are 0.316 it 0.064 
for the bright sample and 0.356 it 0.069 for the faint one, i.e. the widths are the same. A 
Kolmogorov-Smirnov test between the two distributions after shifting all faint bursts down 
in log(r5o) by an amount equal to the difference in the sample means yields a probability of 
0.79 that the two samples are identical. The unchanged width and the high K-S probability 
of equality indicate that the duration difference between faint and bright bursts is achieved 
by an overall shift of the distributions without a change in shape, rather than by a removal 
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of short bursts and a narrowing of the distribution. To further illustrate this, we show in 
Fig. 3 the dependence of the mean and width of the duration distribution on peak count 
rate for all ten groups selected in the above way. The dashed lines are predictions for the 
relations for a standard $7 = 1 universe which has a rate density of gamma-ray bursts that 
is constant per unit comoving volume and time. The photon spectra of gamma-ray bursts 
are assumed to be power laws of slope —1 (Yi & Mao 1994), which matches the average 
observed slope reasonably well. The vertical scale of the curves can be adjusted arbitrarily 
to match the data. The errors are substantial, but the figure does indicate that faint bursts 
are longer without spanning a narrower range of durations, i.e., that the cosmological 
explanation of the time dilation is preferred. 

It should be noted again that since the shifts in both models are small relative to the 
overall width of the duration distribution, it is quite important to account properly for 
selection effects causing incompleteness of the duration distribution, such as the strong 
bias against triggering very short bursts. The test is best done on a sample which has a 
clear maximum in its duration distribution, with the decrease on both sides being real 
rather than caused by selection effects. We believe that the sample we used satisfies 
these requirements, even though the choices we made to define our sample are somewhat 
arbitrary. Small modifications of the sample do not change the result: if we truncate the 
distributions at 1 s, the difference in log(r5o) becomes 0.27 it 0.14. The use of peak count 
rates in 256 ms bins gives the same difference (0.35 it 0.13) as the 64 ms bins, but use of 
1024ms bins or Tgo strongly reduces the difference (typically to 0.12 it 0.12), as expected 
for reasons mentioned above. 

It may be helpful to adopt less arbitrary sample definitions, perhaps by selecting 
on other burst properties. An example would be to exploit the bimodality discovered in 
gamma-ray burst properties (Dezalay et al. 1992, Klebesadel 1992, Kouveliotou et al. 1993) 
and select only the soft group, of which the duration distribution clearly peaks within the 
range of durations that is well-observed by BATSE. The test can no doubt also be improved 
by using the actual BATSE time profiles rather than just the summary data from the 
catalog (see Norris et al. 1994 for a thorough discussion). 

A possible complication that must be kept in mind is that while standard-candle bursts 
are sufficient to model cosmological bursts, a broad luminosity function is also possible. In 
the distribution of peak flux values the effect of that is hardly noticeable in most cases, 
because the distribution is dominated by either the intrisically least luminous bursts or the 
most luminous ones (see, e.g., Ulmer & Wijers 1994). In all these cases, no complications 
in the duration distribution due to intrinsic effects are expected, and our simple picture 
will still be valid. These are also the only cosmological models that matter in the current 
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context, because they are the only ones for which the redshift of the faintest bursts deduced 
from time dilation and that deduced from fitting the peak flux distribution are the same, as 
is the case in the real data. In case the luminosity function is broad, there is a contribution 
from signiflcantly different redshifts at each flux. This would yield extra broadening or 
narrowing of the duration distribution, depending on whether bright bursts are longer or 
shorter than average, respectively. Our test would then fail, but, as stated, we would notice 
this failure because the redshifts of the faintest gamma-ray bursts as deduced from fltting 
standard-candle models to the flux distribution would not be the same as that deduced 
from the time dilation. 

4. Conclusion 

We have presented a method to distinguish the simplest cosmological and local 
interpretations of the observed time dilation of faint gamma-ray bursts relative to bright 
ones. It is based on the fact that the distribution of durations changes with flux in 
very different ways between the two cases. For cosmological bursts, the entire duration 
distribution is shifted to longer durations for weak bursts, and the distribution becomes 
wider by a factor 1 -\- z (hence the distribution of log(r) has a constant width). For local 
bursts, the average duration increases by removing short bursts from the sample while 
leaving the long ones unaffected. This means that the distribution becomes narrower and 
its long-duration end does not shift up in duration. 

The proposed test was applied to a subsample of the second BATSE catalog, and we 
tentatively conclude that the observed time dilation is cosmological in origin. The effect is 
quite signiflcant in our test, and we tried to eliminate selection effects and incompleteness as 
best we could using the summary data in the second BATSE catalog. Given the importance 
of the issue, we nonetheless caution again that the test can be greatly improved by using the 
complete burst time proflles in the manner of Norris et al. (1994) to better deflne samples 
and test selection effects. Also, since many more bursts have now been found by BATSE 
than are included in the second catalog, our test can be applied to an independent or larger 
data set to check that our result is not spurious. 

This work was supported in part by NASA Grant NAG 5-1901. RW is supported 
by a Compton Fellowship (grant GRO/PFP-91-26). This research has made use of data 
obtained through the Compton Gamma Ray Observatory Science Support Center Online 
Service, provided by the NASA-Goddard Space Flight Center. 
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Fig. 1. — Illustration of how the observed log7V(> F) — log arises from a convolution 
of density and luminosity distributions. Diagonal lines in the main panel represent shells 
of constant distance. For distances greater than -Rcore (dashed lines at top left) shells are 
progressively less populated, implying that samples of bursts with fluxes less than F2 have 
relatively fewer intrinsically luminous bursts. This in turn means that the distribution of 
any quantity that is correlated with intrinsic luminosity differs in a sample of bursts that 
are brighter than F2 from a sample of bursts that are fainter than F2. 
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Fig. 2. — Duration comparison of a bright and faint sample of bursts from the second BATSE 
catalog. The brightness measure used is the peak count rate in 64 ms bins; the bright sample 
(solid line) consists of the 10% brightest bursts, and the faint sample (dashed line) consists 
of the bursts that rank between 10% and 20% from the bottom in a brightness-ordered list 
(each sample has 19 bursts). The bin width is log;^Q(2), so a duration change of a factor 
2 means a one-bin shift between the two histograms. The dotted line marks the duration 
cutoff applied to both samples. 
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Fig. 3. — The mean and width of the distribution of log(r5o) for 10 groups of peak count rate 
in 64ms bins as a function of peak count rate. Only bursts with T^q > 3s were used. The 
curves are the predictions (arbitrarily adjusted vertically, and assuming the faintest-but-one 
group has z = 1) for a standard $7 = 1 cosmology with a constant number of gamma-ray 
bursts per unit comoving volume and time. The bursts are assumed to be standard candles 
with photon spectra that are power laws with slope —1. 



